GAUSSIAN TYPE BOUNDS FOR THE NEUMANN-GREEN FUNCTION OF A 

GENERAL PARABOLIC OPERATOR 



MOURAD CHOULLI AND LAURENT KAYSER 

Abstract. Based on the fact that the Neumann-Green function can be constructed as a perturbation of the 
fundamental solution by a single-layer potential, we establish a gaussian lower bound and a gaussian type 
upper bound for the Neumann-Green function for a general parabolic operator. We build our analysis on 
old tools coming from the construction of a fundamental solution of a general parabolic operator by means 
of the so-called parametrix method. At the same time we provide a simple proof for the gaussian two-sided 
bounds for the fundamental solution. 
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1. Introduction 

Let f2 be a bounded domain of K" with '"-boundary, where < a < 1. Let to < ii, set Q = x {to, ti) 
and consider the second order differential operator 

L = aij{x, t)dfj + bk{x, t)dk + c{x, t) ~ dt- 

Here and henceforth wc use the usual Einstein's summation convention. 
We make the following assumptions on the coefficients of L: 

(i) the matrix {aij{x, t)) is symmetric for any [x, t) g Q, 
(ti) a,, g W^^°°{Q), bk, CG C([io,ti],C"(rT)), 
{tii) a,j(x)e,ej > ix,t) e 0, C e M", 

(iv) ||a,y||M/i,o=(Q) + ||&fe||L-(Q) + I|c||l-(q) < A, 
where A > and A > arc two given constants. 

Since we will use the fundamental solution in the whole space, we begin by extending the coefficients of L 
in a neighborhood £7 of O to coefficients having the same regularity. We observe that this is possible in view 
of the regularity of J7. For simplicity we keep the same symbols for the extended coefficients. We may also 
assume that the ellipticity condition holds for the extended coefficients with the same A. Pick tp G C^{n) 
satisfying < < I and t/i = 1 in a neighborhood of fl. We set 

ciij = a^jip + A(5y(l - V'), bk = bkip, c = cip 
1 
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and 

Z = aij{x,t)df^ + bk{x,t)dk +c{x,t) - dt- 

Clearly the coefficients of L satisfy the same assumptions as those of L. So in the sequel we will use the 
same symbol L for L or its extension L. 

We are interested in gaussian two-sided bounds for the Neumann-Green function associated to the operator 
L. More specifically, denoting by G the Neumann-Green function for L, we want to prove an estimate of the 
form 

[C(i-r)]-"/2e-^^ <G(x,<;C,r) < [C{t - T)]-"/^e-^^ , {x,t; ^,t) e {t > r}, 

where the constants C and C depend only on fl, to, ti and A. 

We succeed in proving that the above gaussian lower bound holds true. But we are only able to prove an 
upper bound which is weaker than a gaussian upper bound. Namely we prove an upper bound of the form: 

{{t-Ty/^ + \x-^\)G{x,t,(,T) < [C{t-T)r"/'e-^^, {x,t-^,T)eQ\ t>T. 

for some constant C that can depend only on fl, to, ti and A. 

It is an open problem to know whether the gaussian upper bound is true for a general smooth bounded 
domain Vt. We will see in Section 4 how the two-sided gaussian bounds, for the Neumann-Green function, 
can be obtained in a straightforward way from the gaussian two-sided bounds for the fundamental solution 
when f2 is a half space. 

To our knowledge these kind of gaussian estimates have never been established before. The situation 
is completely different for the Dirichlet- Green function since this later vanishes on the boundary. One 
can prove in a straightforward manner, with the help of the maximum principle, that the Dirichlct-Green 
function is non negative and dominated pointwise by the fundamental solution and so it has a gaussian upper 
bound. Aronson jAr2j (Theorem 8 in page 670 j3 get an interior gaussian lower bound for the Dirichlet-Green 
function. Later Cho |Ch| . Cho, Kim and Park jCKPj extended this result to a global weighted gaussian 
bounds involving the distance to the boundary. 

When L has time-independent coefficients a fundamental solution or a Green function is reduced to a 
heat kernel. We mention that there is a wide literature dealing with gaussian bounds for heat kernels. We 
quote the following classical books: |Da| . |Grj . |Ouj [Sa] . [St], but of course there are many other references 
on the subject. 

As we said in the summary, the main ingredient in our analysis relies on the classical construction of the 
fundamental solution by the so-called parametrix method. We revisit this construction in the next section 
and we derive from it the gaussian two-sided bounds for the fundamental solution. Special attention is paid 
to the dependence of the gaussian upper bound on the lower order coefficients of L. This dependence will 
be in the heart of the proof of the gaussian type upper bound for the Neumann-Green function since we 
will conjugate L with e'^, for some function ij} appropriately chosen. This is a classical argument to derive 
a gaussian upper bound from a Nash upper bound. Note however that this argument is not enough to get 
a gaussian upper bound for the Neumann-Green function. In our case we need to conjugate again with a 
special function together with a maximum principle argument in order to get a gaussian type upper bound. 
This is done in Section 4. Before, we prove in Section 3 a gaussian lower bound for the Ncumann-Grcen 
function. To do so, we construct the Neumann-Green function as a perturbation of the fundamental solution 
by a single-layer potential. The gaussian lower bound is derived from an estimate for the kernel of the 
single-layer potential which is the key point in the proof. 



^Let us observe that Theorem 8 in page 670 of IAr2l can be used to extend the results of Section 3 of IFSI to a general 
parabolic operator. In other words, one can obtain a proof of a continuity theorem by Nash [Na) and the Moser's Harnack 
inequality IMol for a general parabolic operator, which is based on the two-sided gaussian bounds for the fundamental solution. 
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2. The parametrix method revisited 



This section is concerned with gaussian two-sided bounds for the fundamental solution. For a systematic 
study of the fundamental solution we refer to the classical monographs by A. Friedman and O. A. 
Ladyzhenskaja, V. A. Solonnikov and N. N. Ural'tzeva |LSUj . 

In the sequel P = K" x (to,<i). 

We recall that a fundamental solution of Lu = in P is a function E{x, t; ^, r) which is C^'^ in P^r\{t > r}, 
satisfies 

LE{- t) = OinM" x {t < t < ti}, for any r) € R" x [to, hi 
and, for any / e Co(M")0, 



hm / E{x, t- ^, r)/(e)de = fix), x G M". 

The construction of a fundamental solution by the so-called parametrix method was initiated by E. E. 
Levi [Lej . Let a = (a*-') be the inverse matrix of (flij), \a\ the determinant of a and 

/„ /— o(S.T)(x--i;).(x^o 



Z(a;,i;C,T) = [47r(i-T)]-"/V|a|e , (a;,i;^,T) £P'n{t> r}. 

This function is called the parametrix. It satisfies 

LoZ(-,-,^,r) = OinM" X {r < t < h} for any r) e x [to,ii[, 

where 

(2.1) Lo = a,,{i,T)dl~dt. 

In the parametrix method we seek E, a fundamental solution of Lu = in P, of the form 



(2.2) 



E{x,t-^,T)=Z{x,t-^,T) + 



/ / Z{x,t;'q,a)^{7],a;^,T)d7jda, 
Jt Js." 



where $ is to be determined in order to satisfy LE{- ,■ ;(,,t) = for any (^, t) e M" x [to, ti[ 
Following Formulas (4.4) and (4.5) in page 14 of [K], $ is given by the series 



e=i 



where $i(a;, t; ^, r) = LZ{x, t; ^, t) and 



$^+i(x,t;C,T) = /" /" $i(x,t;ry,cr)$f(ry,cr;C, 



T)dr]da, i>l. 



Here for simplicity we write LZ{x, t; ^, r) instead of [LZ{- , •, ^, T)](a;, t) 
Let di, 1 < j < given by 

a'-'"(e,T)(.T,-e,) 



= di{x,t;^,T) = -- 



2(t-r) 



(x,t;e,T)GP^n{t>T}. 



Then 



-^^^ + dd- 
2{t-r) + ^'^^ 



Therefore, taking into account p.ip . we get 

LZ ^ LZ - LqZ = I {aij{x,t) - ay(^,T)) 



"One can take a larger class of functions. Namely a class of continuous functions satisfying a certain growth condition at 
infinity (e.g. formulas (6.1) and (6.2) in page 22 of IPrp. 
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We write LZ = '^Z, where 



+ bkdk + c. 



Let 



M = max||ay||^i,^(Q), N = max(max ||6a;||loc(q), ||c||ioo(Q), 1). 

i.j k 



Since 



we have 
(2.3) 



\d,\ < 



M\x-^\ 



t-T ' 

\a,j{x, t) - a,y(C, t)| < A/(|.T - el + t - r), 



Here P is a polynomial function of degree less or equal to three whose coefficients depend only on M . 
Unless otherwise stated, all the constants we use now do not depend on TV. 
In light of (|2.3p we obtain 



\LZ\ < CN{t-T) 



-(n+l)/2 



P{Tj)e 



with 



-(A/4)|r,|2 



-(A/8)|r,|2 



But the function p £ (0,+oo) — > P{p)e^''^/^^P is bounded. Consequently, where A* = A/8, 



(2.4) 



|a>i(x, r)| = |LZ(z, r)| < iVC(t - r)-("+i)/2e-^ 



The following lemma will be useful in the sequel. Its proof is given in page 15 of [F?] 
Lemma 2.1. Let c > and — oo < j,/] < n/2 + 1. Then 



{t — a) '^e t-" {a - r) '^e 



' drjda 



/a \ ^/'^ 

i JL \ B{nl2 - 7 + 1, n/2 - /? + l)(i - r)"/2+i-7-/3e-- 



where B is the usual beta function. 
We want to show 



(2.5) 



\^i{x,t;^,T)\ < {NCY&-\t-r)-^^+'-'y'l[B{l/2,j/2)e- 

Here C is the same constant as in ()2.4p and C ~ (fr)"^^- 
As 



> 2. 



<^2{x,t;tT) ^ [ I <S>i{x,t;r],a)<i>i{T],a;^,T)dijda, 



estimate (|2.4[) and Lemma [2TT] with j — /3 = n/2 + 1 show that (|2.5p holds true with £ — 2. The general case 
follows by an induction argument in £. Indeed, using 

/■* 

<S>i+i{x,t;£„T) 



/ / '^i{x,t;ri,a)^i{r],a;^,T)d'qda, 



(123), for f and Lemma [Q with 7 = n/2 + 1 and /3 = (n + 2 - £)/2, we obtain easily that (^3]) holds 
true with ^ + 1 in place of ^. 
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If r is the usual gamma function, wc recall that 

S(1/2,j72) 



e-i 



i>i 



r(i/2)r(j72) 
r((j + i)/2)- 

T{1/2Y V 



T{£/2) T{£/2)' 



Therefore 

(2.6) Y[b{1/2,j/2) 

([2^ - ([2^ entail that 

(2.7) |<i>(x,t;^,r)| <^|$,(x,t;e,r)| <C(l + 5)(t-r)-("+i 

with 



5 = ^ [cA^(t-T)i/2j /r{e/2). 



e>i 



We have r(£/2) = r(m) = (m - 1)! if £ = 2m and r(^/2) = r(7n + 1/2) > r(TO) = (?7z - 1)! if ^ = 2m + 1. 
Then 



_ — _ r -| 2m _ — _ r 



2m+l 



7r(2m + l/2) 



m>0 



2m+l 



/(m- 1)!. 



Whence 



Plugging this estimate into p.7p . we obtain 



(2.8) 



|$(a:,<;e,T)| <C(t-T) 



-(n+l)/2. 



-+CAr2(t-r) 



With the help of Lemma 12.11 estimate ()2.8p yields 



(2.9) 



Z(x, i; 77, ct)<I>(77, cr; ^, T)dr/da 



< 



We obtain as a consequence of this inequality 

(2.10) \E{x, t- t)| < C{t ~ r)-"/2e-^%^+CAr=(t-r)^ 

This estimate is essential when establishing the gaussian type upper bound for the Neumann-Green function. 

In the rest of this section wc forsake the explicit dependence on A^. So the constants below may depend 
on f7, A, A, Iq and ti. 

From (|2.9p we deduce in a straightforward manner 

(2.11) E{x,t;^,T)>d{t^T)-"/^, {x,t;tT) e P^, t> T, C\x-^\^ <t-T. 

By Theorem 11 in page 44 of [E], E is positive. Moreover E satisfies the following identity, usually called 
the reproducing property, 

(2.12) E{x,t;^,T)= E{x,t;ri,a)E{T],a;^,T)d7j, x,CeR", to<T<a<t<ti. 

We can now paraphrase the proof of Theorem 2.7 in page 334 of [FS| to get the gaussian lower bound for E. 
We sum up our analysis in the following theorem. 

Theorem 2.1. The fundamental solution E satisfies the gaussian two-sided bounds: 

(2.13) [C(t-T)]-"/2e-c^ <E{x,t;^,T) < [C{t ~ T)]-''/^e-^^ , {x,t,^,T) e x {t > r}. 
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Remark 2.1. Let us observe that the proof of Theorem 12.11 presented here is elementary in comparison with 
that given in [FS| for proving the gaussian two-sided bounds for the operator di{aij{x,t)dj ■)—dt with (C°°) 
smooth coefficients. Of course, Theorem 12.11 applies the operators in divergence form. Note however that 
the constants C and C appearing in p.l3|) do not depend on to and ti in the gaussian two-sided bounds 
established in |FSj . We mention that gaussian two-sided bounds were obtained by S. D. Eidel'man and F. 
O. Porper [EP| when the coefficients of L satisfy the uniform Dini condition with respect to x. The main 
tool in |EP| is a parabolic Harnack inequality. We refer also to jArlj . [Fa] , p] and [NSj . where the reader 
can find various results on bounds for the fundamental solution. 

3. Gaussian lower bound for the Neumann-Green function 

The unit outward normal vector at a: G is denoted hy v = i'{x). Henceforth S = 951 x [tQ^ti). 

For T £ [to,ti[, set ~ fl x (t, ti) and Y.^ = x (r, fi). We consider the Neumann initial-boundary 
value problem (IB VP in short) for the operator L: 



From Theorem 2 in page 144 of and its proof, for any ip G C^{^1), the IB VP (|3.ip has a unique 
solution u G C{Qt) n C^'^(Qr) given by 



(3.1) 




(3.2) 




Here 



(3.3) 




with 




Ml = 2d^E 




Let, where {x,t) e St and £, G il, 




Assume for the moment (see the proof below) that 



(3.4) 




We set 



(3.5) 




)dr/da + E{x, t] ^, r). 



It follows from Fubini's theorem that 



(3.6) 




The function G is called the Neumann-Green function for Lu = in Q. 
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Wc have, for any < ip € C^{Q), u > according to the maximum principle (e.g. Theorem 2.9 and 
remarks following it in page 15 of [Li]) ; whence G > 0. From the uniqueness of the solution of the IBVP 
(13.11) wc have also 

G{x,t;^,T)'ilj{£,)d^ ^ [ G{x,t,7],a)d7] [ G{r,,cj,^,T)i^{£)d^ for any i; e C^^ {Q), t < a < t. 
in Jn Jn 

Therefore 

(3.7) G{x, t; r) = / G{x, t, r], a)G{ij, a, ^ T)dij, t < a < t. 



That is G has the reproducing property. 

Let us observe that, when c — 0, G satisfies in addition 



/ Gix,t;^,T)d^^l. 
Jn 



We shall need the following key lemma for establishing the gaussian lower bound for G. 
Lemma 3.1. For 1 — a/2 < ^ < 1, we have 

(3.8) |7V(x,i;C,T)| <C(t-T)-'^|x-^|-"+\ (x, i) € E„ ^ € f}, x ^ ^. 

The lemma below appears in page 137 of [ft] as Lemma 1. It is needed for proving Lemma |3. II 
Lemma 3.2. Let < a,b < n — 1. Then 

(3.9) / \x-T]\-''\7j-^\-^df]<C\x-^\''-^-^''+''l 

Jdn 

Proof of Lemma WJ\ Let 1 — a/2 < /i < 1 be given. From formula (2.12) in page 137 of |ft], we have 

\d,E{x,t-S„T)\ < C(i-r)-^|a;-er"+'+^"'+""'^ 

and then 

(3.10) \Mi{x,t;^,T)\ < C(t-T)-'^|a:-^|-"+i+(2'^+"-2). 
Since 



|M2(x,i;C,r)| < / / \Miix,t;f^,a)\\Mi{f^,a;tT)\dvda, 
Jt Jdn 



(Pl^ leads 



(3.11) |M2(x,i;C,r)| < CM {t - a)-^{T - a)-^da |x - r;r"+^+^^'^+"-^^ IC - ??r"+'+^'^+"-''dr;. 

Jt Jdn 

Or from Lemma 13.21 

(3.12) / \X- ^|-"+1 + (2a<+o-2)|^ „ ^|-„+l+(2M+a-2)^^ < _ ^|-„+l+2(2M+a-2)^ 

Jon 

On the other hand 

(3.13) J^{t- ayir - ayda = {t ~ ^ ^-^(i „ sf'^ds ^ {t - r)-^+(i-^)S(l - /i, 1 - /^)- 

We plug (jXT^ and (|XT^ into ([XTT]) ; we obtain 

|M2(x,t;C,T)| < C2C(t-T)-'^+(i-^)B(l -^,1 -/i)|a;-^r"+i+2(2A'+"-2). 

Now an induction argument in £ yields 

\Me{x,t;tr)\ < C'C'-\t - ^(1 - m)' , _ ^|-„+i+£(2p+a-2)_ 

i (t(l /ijj 
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It follows from T{£(1 — /i)) > [^(1 — /i)]! (here [•] is the entire part) that the series 

converges uniformly in all its arguments. We complete the proof by noting that 

iV(x, t; e, r) = (t - r)-^|x - ^\-"+'S{x, t; r). 



Proof of (1311). Let 

^ pi p 

^ ft f 

= 2F,(x,t)+2E / / Mt{x,t-i,a)Fr{£.,(T)d^d(j. 
In light of Lemma 13.21 and with the help Lebesgue's dominated convergence theorem we can assert that 



□ 



Or according to Funini's theorem 

^k{x,t) = [ Nk{x,t;^,T)^{OdC 
Jn 

But ipk{x^t) — > ip{x,t) when k tends to infinity. Then the uniqueness of the limit yields 

ip{x,t) ^ I Nix,t;^,T)^(Od^. 



□ 



We are now ready to prove 
Theorem 3.1. The Green function G satisfies the gaussian lower bound: 

(3.14) [C{t-T)]-"/^e-^^-^ <G{x,t;^,T), {x,t;^,T) e n {t > t}. 

Proof. Let 



Goix,t;^,T) ^ / E{x,t;Tj,a)N{r],a]^,T)d^da. 
Jt Jdn 

Let < /3 < 1 be fixed. From the gaussian upper bound for E we obtain in a straightforward way that 

|i?(x,t;e,r)| <C(t-r)-^|x-Cr"+2^. 
This estimate together with (j3.8p imply 

|Go(x,i;^,r)| <C f{t-cj)-^\a-T)-"da f \x - v^+'lv - C^+'P di^ 

< C{t-T)^~''-^\x~£_\-"+^'^, 
where we used Lemma [321 That is 

(3.15) |Go(a;,i;e,r)| <C(t-r)-"/2+i-"[(t-r)-i|a;-eP]-"/'+^. 

But we already know from ()2.1ip that 

E{x, t- t) > G{t - t)-"/^ (a;, i; f , r) e P^ t > r, d\x -£.\^ <t-T. 
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Hence 



G(x,t;e,r) > E{x,t;^,T) - \Go{x,t; ^,t)\ 

> C{t - t)-"/2(^ _ (j^j^ _ rf-t^)^ t > r, C\x - < i ■ 



r. 



Consequently, 

G{x, t- r) > C(t - t)-"/2, (a:, <; ^, r) e g^ t > r, C\x -i\^<t-T. 

Or cquivalcntly 

(3.16) G(a;,t;^,T) > C(i-r)-"/2, (x, t; ^, r) e Q^ i > r, |a; - ^1 < - ^)^/^- 

Since G > and satisfies p.7p . we proceed as in the proof of Theorem l2.1l to get the gaussian lower bound for 
G. Here we rewrite the arguments for the reader convenience. As $7 is connected, we find a path 7 : [0, 1] ^ fi 
connecting x to ^ which is piecewise constant. Let fc be a positive integer and set yi = "f{i/k), < i < k. 
Then it is not difficult to show that there exists a constant c > 1 not depending on k such that 

(3.17) \y^+l - y^\ < - 0<i<k-l. 

Shortening if necessary the constant G in p.l6|) , we make the assumption 

(3.18) C{t - t)1/2 < 4dist(7([0, l]),dn). 

When 2c\x - ^| < d(t - t)^/^ (implying |a; - ^| < C(t - t)^^"^), (j3ll)) follows immediately from (l3l6l) . 
Therefore we may assume that 2c\x — ^\ > C{t ~ t)^/^. We choose m > 2 as the smallest integer satisfying 

2ct^<a(t-r)i/2. 

771 ' 

Set Xi = ^{i/m), < i < and 



In light of the reproducing property and the positivity of G we obtairH 

G(x,t;e,r) = ^^...|^G(x,t;a /"'~^^^^^ )...G( ^^^7n~^^^ Cn.-i;^,r^rfei...rfgn.-i 

(3.19) > / G (x, t; -'^' + ^ )...g( '+ - , A d^r ■ ■ ■ d^-^- 



JB{xi,r) JB[xrr,-i,r) 

Let ^0 = 2; and = ^ ; we have 



Whence 



\x - f\ \x - f\ 

- < - Xi\ +2r< cJ ^ + 2r < c ' + 2r < 4r, < i < m - 1. 



'Note that B(xi,r) C 1 < i < m — 1, as a consequence of 113.1811 
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It follows from (|3.16p that 

G{x,t;^,T)> ... -) d^i...dU-i 



, \ — nm/2 



ft 

16 



m 

-| n(m-l)/2 



^/2 



where Wn is the measure of the unit ball of K" . In particular 

(3.20) G(a;, <; r) > Ce"^™(< - t)-i/^ 
Or from the definition of m we have 

(3.21) ™-i<yL_^. 

Finally, a combination of (|3^ and ((3?2T|) leads to (|3TT4)) . □ 

We mentioned previously that, when i has time-independent coefficients, the Neumann-Green function in 
nothing else but the heat kernel associated to the semi-group generated by L under the Neumann boundary 
condition. 

To this end we assume that is of class Q. to = 0, ii = T and the coefficients of the operator L are 
time- independent. In other words L is of the form 

(3.22) L = a,,{x)dl + bk{x)dk + c{x) - dt 
and the following assumption holds true: 

[i') the matrix (aij{x)) is symmetric for any x £0,, 
(^^') € W^'^{Vl), bk, cG C"(f7), 
(in') a.,(x)60 > {x,t) e n, C e M", 

{iv') ||a,y||n/i,oo(o) + WbkWL'^in) + \\c\\L'^(n) < A, 
where A > and A > are two given constants. 

Under the above mentioned assumptions, it is known (e.g. |0u| ) that the operator 

Al = aij{x)dfj + bk{x)dk + c(.t), with domain D{Al) = {u e H^{il); dyU = on dVL}^ 

generates an analytic semi-group e*'^^ on i^(fi). 

Let -0 G C^iP)- Since u{t) = 6*^^-0 is the solution of the IBVP ((3T1) . 



6*^^^^= / G(a;,t;C,0)V'(O^^C, 0<i<r. 

We rewrite this equality as follows: 

e*^^V= / K{x,£„t)iij{£,)di, 0<t<T. 
n 



^This assumption is not really necessary, we make it just because in this case the domain of the operator Aj^ below is a 
subset of H^{Q,) and therefore the normal derivative of an element of D{Ai^) exists in the usual trace sense. 
^Notice that A]^ can be rewritten as an operator in divergence form simply by observing that 

aij{x)dfj = dj{aij{x)dj- ) - djaij{x)dj. 
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The function 

is usually called the heat kernel of the semi-group e*^'^ . 
A straightforward consequence of Theorems 13.11 is 

Corollary 3.1. The Neumann heat kernel K satisfies the gaussian lower bound: 
(3.23) (Ci)""/^e-^^^ < K{x, ^, t), (x, € < t < T. 

The gaussian lower bound for the Neumann heat kernel was proved in |COY] when L is the laplacien. A 
quick examination of the proof show that this result can be extended to an operator of the form p.22p with 
C°°-smooth coefficients. The key point is the Holder continuity of x — > K{x, ^, t) which relies on the fact 
that fi — A L is a.n isomorphism from H^{il.) into H^~'^ for some large /i and s of order n/2 + 1. This 
explain why the approach in [COY] can not be used to extend the lower gaussian bound for operators with 
less smooth coefficients. 

4. Gaussian type upper bound for the Neumann-Green function 

In this section, il is assumed to be of class C^'". 
We start with the so-called Nash upper bound. 

Lemma 4.1. The Green function G satisfies the Nash upper bound: 

(4.1) G(x,t;e,T) <C(t-T)-"/2, {x,t;^,T)eQ',t>T. 
Proof. From the preceding proof we have, for any < /3 < 1, 

(4.2) G{x, t; e, r) < G{t - r)-^|.T - ^\-"+^^ . 
Or by (|3.7p (reproducing property) 

(4.3) G(x, t; v,t)= f G{x, t, C, {t - T)/2)G{t {t - t)/2, T)d^. 



In light of the inequality in Lemma 2 in page 14 of jFr)^ . it follows from (|4.2p and (|4.3p 

G(x,t;C,r)<C(t-r)-^^|x-^r"+2(2«. 
We find by repeating this argument fc-times 

(4.4) G{x,t;^,T) < G(i-T)-'='5|a;-^r"+*^(2/3)^ 

Let k be the smallest integer satisfying n < 2k. Then the choice /? = n/{2k) in (14.41) gives (j4.ip . □ 

Let ( e R" and set |i/|oo = HIj^HIoo- As il is of class C^ ", we obtain from Lemma 3.1 in |Choj that there 
exists 6*^ e C'^'"(n) having the properties 

(4.5) o^>i, e( = iondn, -d^e^: > \iy\oc\C\ on on. 

Moreover examining the proof of Lemma 3.1 in |Choj we obtain in a straightforward manner that 

(4.6) \W, |%0cl < e\a 

where p is a constant independent on 

We now conjugate L by e^'^+'^'^l *6'^, where K is a. constant to be specified later. To this end we consider 



''Let < a, (3 < n. Then there exists a positive constant C such that for any (x,^) g fP, a; ^ we have 

/ \x ~ ri\~"\rj~ ^l'/^ dri < C\x ~ if a + 13 > n and f |x - ryp" |r? - ^p'^dr? < C if a -I- /3 < n. 

Jn Jn 
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An elementary calculation gives L« = L« + c« with 

- a,,dl + bpdp + 2akiiCi + B-^diB(^)dk + 6„,(C,„ + B-^d„,ec) + c + 9^'a,,d^^9c - dt 

and 

In light of (|4.5p and (|4.6p we can choose A', independent on C, such that c*" < for any ^ e R". 

Let < V e C^{^) and li the corresponding solution of the IBVP Then v ^ e-'^-^--^^l^l'*6'~^w 

solves the following IBVP 

r L^w = in Qr, 

(4.7) i w(- ,t) = e-^-^--^lf'^6l^Vj inf^, 

Next, let V be the solution of the IBVP 

r Z'^u = iuQ^, 

(4.8) I v(-,t) =e-«-^-^lfl'^6'^V infl, 

I. d^v = on S,- . 

Since v > 0, 

L?(w- w) = -c^-y > iuQ^, 
(i;-w)(-,t) =0 inrj, 
^^(w- w) = -(C- i' + (9^6'c)w > onE^. 

The last inequality is obtained from the third inequality in (|4.5p . We invoke Theorem 2.9 [Li] to deduce that 

V < V. 

Let E'' and be respectively the fundamental solution and the Neumann-Green function corresponding 
to L''. It follows from estimate p.lOp that 

EHx,t;^,T) < C(i-r)-"/2e-^4^+'5(i+ICI+ICn(t-r), 
A slight modification of the proofs of Lemmas 13.11 and 14.11 leads to the following Nash upper bound for G"" : 

(4.9) G';(a:,t;e,T) <C(t-r)-"/^e'5(i+l«l+lfl')(*-^\ 
Therefore, for any < V e G(f (17), 

Or 

uix,t)= f G{x,t;^,T)i;iOdC 
Jn 

Whence 

Gix,t,^,T) < e«-+^l«l'*G?(x,t;^,r)e-«-^l?l'-0c(x)0^-i(C). 

This and ((49| imply 

(4.10) G(a;,t,C,T) < G(i - r)-"/2gC-(:'=-e)+c(i+|C|+|CI^)(t-r)g,^(^.)g,-i 
and then 

G(x,t,e,T) < G|C|(t - r)-"/2e«-(^-«)+^(i+l«l+l^l')(*-^). 
Taking — x)/[2C{t — r)] in the last estimate we get the following result. 

Theorem 4.1. The Neumann- Green function satisfies the upper bound: 

((i-r)i/2 + |2;-ei)G(x,i,C,r) <[G(t-r)]-"/2e-^fe^, {x,t; ^,t) e t > r. 
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Remark 4.1. 1) Returning to the proof of Lemma 3.1 in [Choj . we prove that for any compact subset K of 
n there exists Ck such that 

\0cix)e^HO\ <Ck, xeTl, (,eK and c e M". 

We derive from (|4.10p the following gaussian upper bound 

G(x, t, r) < - T)-"/2e-^^fe^ , {x, t; ^,t) e n e K, t > r}. 

2) The construction of the Neumann- Green function is quite simple when is a half space, say 

= M!^ = {x = (.t',.t„) e M""i X M; Xn > 0}. 

Assume for simplicity that the original coefficients of L satisfy assumptions (i)-(iv) in the whole space R" 
and set 

G{x, t; tT) = - [E{x, t; r) + E{x',t; ^, r)] , x^^eM.^, T,te {to,t,), t < t, 

where = (x' , — a;„) if a; = (a;', a;„). 

It is not hard to check that G is the Neumann-Green function for L in M" x (tg, ti). In addition observing 
that |x — ^1 < jx" — ^1 for any a;,^ € M", we see the gaussian two-sided bounds for E yield the gaussian 
two-sided bounds for G. 

For 1 < p,q< oo, let j| • \\p,q denotes the norm in L'^{to,ti; LP{^l)). We have similarly to Corollary 7.1 in 
page 668 of jAr2j the following result which is a direct consequence of the upper bound in Theorem 14.11 

Corollary 4.1. Let 1 < p,q < oo be such that 

n 2 

- + - < 1. 

p q 

Then 

\\Gix,t;-,-)\\p,^,,, ||G(•,•;^,T)||p^,- <C, 
where p',q' are the respective Holder conjugate exponents of p and q. 

Let / e C^{Qr) and p, q be as in the previous lemma. Then the solution of the following BVP 

Lu = / in Qr, 
u{- ,t) ~ in 17, 
d^u = on El- 
can be represented by the formula 

uix,t)= f f G{x,t;^,T)fi^,T)dCdT. 
Jt Jn 

In light of this formula we have the following immediate consequence of Corollary 14.11 

||lt||oo,oo < CWfWp^q. 
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